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Приведены результаты исследований, позволяющие установить влияние структурных особенностей образцов на характер вол- 
новых процессов в них. С позиций механики сплошных сред проводится моделирование волновых процессов при действии им- 
пульсной нагрузки на часть поверхности образца. Численный анализ показывает связь характеристик неоднородностей с изме- 
нением параметров волновых процессов, протекающих в твердотельной среде. 


Для проведения неразрушающего контроля 
композиционных материалов и мониторинга изме- 
нений напряженно-деформированного состояния 
(НДС) массивов горных пород в шахтных условиях 
разрабатывается метод, использующий эффект ме- 
ханоэлектрических преобразований в диэлектри- 
ческих материалах [1,2]. Схематически эти преоб- 
разования можно объяснить следующим образом. 
В композиционных материалах всегда присутству- 
ют структурные неоднородности и границы разде- 
ла фаз, составляющих образцы. В материалах есте- 
ственного происхождения неоднородности связа- 
ны с наличием включений из других минералов и 


пород или дефектов в виде трещин или пустот. При 
действии изменяющейся нагрузки материал де- 
формируется, при этом границы раздела фаз или 
разного рода неоднородностей являются зонами, 
где наиболее резко выражены градиенты переме- 
щений, деформаций и напряжений [3, 4]. Кроме 
того, на этих границах возникают двойные элек- 
трические слои, образование которых обусловлено 
накоплением и разделением зарядов разного знака. 
Движение зарядов или изменение дипольных мо- 
ментов этих слоев в поле акустических колебаний, 
возникающих при деформировании материала, 
приводит к появлению электромагнитного откли- 
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ка. В результате использования внешнего детерми- 
нированного акустического воздействия на гетеро- 
генные материалы возможно проведение неразру- 
шающей диагностики структуры материала и изме- 
нения НДС массивов горных пород [2]. 

Возникающие при исследовании этого явления 
задачи можно разделить на прямые и обратные. В 
первом случае речь идет о сопоставлении электро- 
магнитного отклика нагружаемого образца или из- 
делия с известной структурой и заданного внешне- 
го силового воздействия. Во втором случае по элек- 
тромагнитному отклику на заданное силовое воз- 
действие необходимо определить размеры, форму 
и характер неоднородности в материале образца 
или изделия - собственно, это и является задачей 
дефектоскопии. В настоящей работе рассматрива- 
ются прямые задачи, которые относительно про- 
сты, но их решение является необходимым этапом 
решения задач второго типа, более сложных и бо- 
лее важных в прикладном отношении. 

Целью работы являлось численное моделирова- 
ние и исследование с использованием аппарата ме- 
ханики сплошных сред процессов распростране- 
ния упругих волн в неоднородном образце при им- 
пульсном воздействии на статически нагруженные 
лабораторные образцы и закрепленные различным 
способом. 

Результаты такого моделирования будут полезны 
при анализе параметров электромагнитного отклика 
статически нагруженных диэлектрических образцов 
при изменении характеристик импульсной нагрузки, 
геометрии и свойств включений. Они помогут уста- 
новить связь параметров электромагнитной эмиссии 
с изменениями напряженно-деформированного со- 
стояния образцов и массивов горных пород. 

В лабораторных исследованиях использовались 
образцы размером 50x50x100 мм 3 . Статическое од- 
ноосное нагружение осуществлялось вдоль наи- 
большего размера. Материалом служил цементный 
камень, который являлся однородным по составу 


диэлектриком. Неоднородности в лабораторных 
условиях создавались в виде включений, для кото- 
рых известны форма, размеры и свойства. Включе- 
ния, имеющие вид параллелепипеда, создавали из 
дерева и меди, при этом грани вставки ориентиро- 
вались параллельно граням образца. 

Эксперимент реализовывался в двух вариантах 
(рис. 1). В первом случае верхняя и нижняя грани 
образца фиксировались без сколько-нибудь значи- 
тельного поджатия (рис. 1, а). Во втором случае об- 
разец помещался между плитами пресса, которые 
создавали статическую нагрузку ^=соп8І вдоль на- 
ибольшей плоскости образца (рис. 1, Ъ). Это имеет 
целью имитировать начальные напряжения, всегда 
присутствующие в строительных конструкциях 
или горных породах. В обоих случаях импульсная 
нагрузка прикладывались к части боковой поверх- 
ности одной из вертикальных граней образца. 

Из динамической теории упругости известно, 
что для однородных сред неограниченной протя- 
женности анализ двух типов волн (продольных и 
поперечных) в двух измерениях дает такие же ре- 
зультаты, что и в трехмерном измерении. Поэтому 
для наших экспериментов задача ставилась и чи- 
сленно реализовывалась как двумерная на случай 
плоского деформированного состояния. Этот под- 
ход существенно снижал вычислительные сложно- 
сти, но позволял получать результаты, физически 
оправданные и качественно отражающие реальные 
трехмерные процессы. 

Применительно к рассматриваемому классу за- 
дач наибольший интерес представляли области, где 
параметры НДС меняются наиболее интенсивно, 
т.е. градиенты соответствующих характеристик 
максимальны. Эти градиенты наиболее выражены 
в окрестностях разного рода неоднородностей, где 
деформационно-прочностные параметры среды 
резко меняются. Решение предполагалось прово- 
дить как в упругой постановке, так и с учетом эф- 
фектов пластичности. 



Рис. 1 . Схемы приложения нагрузки: а) неподвижная фиксация образца; Ь) предварительное поджатие образца с включени- 
ем; сила приложена на боковой поверхности 
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С позиций механики деформируемого твердого 
тела упругая задача формулируется на случай пло- 
ского деформированного состояния следующим 
образом. В декартовых координатах х, у уравнения 
движения имеют вид 


Р 



да дет 

хх _| л У_ 

дх ду 


д 2 ѵ да да 

Р Т = — + — 

ді 2 дх ду 


Здесь р=р(х,у) - плотность материала, 
и=и(х,у,і), ѵ=ѵ(х,у,і) - смещения точек среды 
вдоль осей хи у соответственно, і — время, 
а и =а хх (х,ур), а х =а ху (х,ур), а у =а уу (х,ур) - компо- 
ненты тензора напряжений. Деформации связаны 
с ними соотношениями [5]: 



а 


XX 


-у, а 

1 УУ 



ч 




-у, а 

1 XX 



? 




Здесь V! = , Е { = г , С = где, 

1 — ѵ 1 — ѵ 2(1 + ѵ) 

в свою очередь, ѵ - коэффициент Пуассона, 
Е— модуль Юнга, С - модуль сдвига. Наконец, де- 
формации выражаются через перемещения линей- 
ными соотношениями 

ди ду ди дѵ 

в = — , в = — , 8 = \ . 

ох ду ' ду дх 


Краевая задача замыкается постановкой на- 
чальных и граничных условий. 

В качестве начальных принимали условия нена- 
пряженного и недеформированного состояния. В 
этом случае нижний торец образца неподвижен 
у=0 и для него соответствуют условия отсутствия 
смещений и(х,0,/) = ѵ(х,0,/). Для верхнего торца у=1 
(/ - высота образца, рис. 1) смещения в горизон- 
тальном направлении также отсутствуют и 
и(х,1,і)= 0. При первом варианте проведения экспе- 
римента (рис. 1, а), когда образец просто зафикси- 
рован в осевом направлении, смещения вдоль оси у 
отсутствуют и ѵ(х,/Д)= 0, а во втором варианте экс- 
перимента на верхней грани образца задана поджи- 
мающая нагрузка а уу (х,1,і)=-Р. Правая кромка об- 
разца Х—а свободна от внешних нагрузок 

а хх ( а ,УР) = о ху { а ,у,1) =0 - 


На левой границе х=0 вдоль нормали к ней при- 
ложена импульсная нагрузка, для которой распре- 
деление по поверхности образца в плоском случае 
У=Ѵ(у) по линии х=0 задается колоколообразной 
функцией [6] 


У (у) = ехр 


(у-уоУ 



Изменение импульсной нагрузки во времени 
Т=Т(і) определяется затухающей синусоидой 


т = е 


( 2 пі\ 


81П 


До у 


В результате на границе образца нормальная со- 
ставляющая напряжений определяется зависимостью 


сУуу{0,ур)=У{у) Т {і). 

Касательной составляющей нагрузки на этой 
границе пренебрегаем, т.е. ст уу (0,у,і)=0 для всех зна- 
чений 0 <у<1 

Численная реализация проводилась с использо- 
ванием нецентральной разностной схемы второго 
порядка точности относительно шагов по про- 
странству и времени [7]. Корректность численных 
результатов оценивалась как по внутренней сходи- 
мости результатов при изменении параметров ко- 
нечно-разностной сетки и шагов интегрирования 
по времени, так и проведенными авторами расче- 
тами модельных задач [8, 9]. 

Численные расчеты проведены для образцов 
двух типов. Первый из них представлял собой од- 
нородный образец из цементного камня. Во втором 
случае рассматривались образцы, в которых есть 
включение в форме прямоугольного параллелепи- 
педа. Вдоль оси статического нагружения его раз- 
меры в плане совпадают с поперечными размерами 
всего образца. Таким образом, неоднородность в 
образце представляла собой поперечный слой с ха- 
рактеристиками, отличными от основного матери- 
ала. Импульсная нагрузка прикладывалась либо к 
ненапряженному образцу, либо к образцу, предва- 
рительно нагруженному с усилием Р. Ниже рассмо- 
трены четыре варианта численных расчетов. 

На рис. 2-4 приведены результаты расчетов в 
виде распределения изолиний интенсивностей де- 
формаций по расчетной области для четырех мо- 
ментов времени і\ а) 5; Ь) 10; с) 15; б) 20 мкс. Тем са- 
мым принимается, что эволюция интенсивности 
деформаций, являющейся своего рода интеграль- 
ной характеристикой деформации образца в каж- 
дой его точке, дает представление о развитии вол- 
нового процесса. Вдоль осей размеры расчетной 
области, согласующиеся с размерами лабораторно- 
го образца, указаны в мм. Шаг между значениями 
интенсивности деформаций, отвечающими сосед- 
ним изолиниям, равен 6Т0 -7 . 

На рис. 2 приведены результаты для однородно- 
го образца, в начальном состоянии не подвержен- 
ного действию внешней нагрузки. В пространстве 
центр ударного импульса расположен от нижнего 
основания на расстоянии 0,75 высоты образца. На 
рис. 2, а , (/=5 мкс) распространяемая зона упругих 
возмущений имеет почти сферически симметрич- 
ную форму. После выхода волны на верхнюю 
(рис. 2, Ь , /=10 мкс) и правую (рис. 2, с, /=15 мкс) 
кромки образца на распространяющуюся волну на- 
кладываются отраженные от этих кромок возмуще- 
ния. К моменту времени /=20 мкс (рис. 2, д) упру- 
гое возмущение распространяется по всей расчет- 
ной области, при этом какие-либо особые зоны в 
образце не просматриваются. 

Приложение предварительной сжимающей на- 
грузки практически не влияет на процесс распростра- 
нения упругих возмущений по образцу. Основные чер- 
ты процесса (рис. 3) почти полностью повторяются. 
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Рис. 2. Поля интенсивности деформаций в однородном образце, в начальном состоянии не подверженного действию внешней 
нагрузки, для последовательных моментов времени 
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Поля интенсивности деформаций в однородном образце, в начальном состоянии подверженного действию внешней 
нагрузки, для последовательных моментов времени 
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Рис. 4. Результаты для неоднородного образца, в начальном состоянии не подверженного действию внешней нагрузки, для 
последовательных моментов времени 
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Рис. 5. 
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Результаты для неоднородного образца, в начальном состоянии подверженного действию внешней нагрузки, для по- 
следовательных моментов времени 


Для тех же моментов времени результаты расче- 
тов приведены на рис. 4, когда в образце имеется 
неоднородность в виде относительно жесткой 
вставки высотой 1 см на расстоянии 5 см от нижне- 
го основания образца. Уже к моменту времени 5 мкс 
упругая волна достигает этой вставки, и начинается 
процесс частичного отражения волны и частичного 
проникания ее внутрь вставки. В последующие мо- 
менты времени по характеру распределения изоли- 
ний деформаций место и размеры положения 
вставки легко отслеживаются. Границы вставки яв- 
ляются зонами, где градиенты деформаций много 
больше, чем в остальной расчетной области. 

Такой же характер распространения и концен- 
трации деформаций прослеживается в целом и при 
моделировании предварительно поджатого в вер- 
тикальном направлении образца (рис. 5). 

В целом анализ результатов показывает, что на- 
личие вставки приводит к возникновению зон кон- 
центрации деформаций и напряжений. Расчеты по- 
казали, что распределения изолиний интенсивно- 
стей напряжений практически совпадают с приве- 
денными результатами для интенсивностей дефор- 
маций. Именно в окрестностях таких границ возни- 
кают наибольшие градиенты параметров НДС и бу- 
дет возникать максимальный поляризационный за- 
ряд. В результате из зоны с повышенным градиен- 
том параметров НДС будет регистрироваться мак- 
симальный по амплитуде электромагнитный сиг- 
нал, возбуждаемый внешним акустическим импуль- 
сом или акустическими сигналами, возникающими 
при образовании трещин в напряженной зоне. 

В ранние моменты времени, когда упругая вол- 
на проходит по образцу, и ее фронт не выходит ни 
на включение, ни на границу расчетной области, 
форма фронта близка к сферической (в плоском 
случае к полуокружности). Это вполне согласуется с 
физическими представлениями о моделируемом 
процессе и подтверждает корректность как самой 


модели, так и метода ее численной реализации. По- 
сле выхода фронта волны на границу раздела основ- 
ного материала и включения или на границу рас- 
четной области проявляются эффекты, связанные с 
отражениями волн, и на приведенных иллюстра- 
циях они видны как сгущения изолиний, отражаю- 
щие большие градиенты соответствующих величин. 

По результатам численного моделирования 
можно сделать следующие выводы: 

1 . Предварительное поджатие образца в виде прило- 
женного одноосного сжимающего напряжения 
практически не сказывается на изменении про- 
цессов распространения упругих волн. Некоторое 
увеличение скорости фронта волны по поджатому 
образцу, которое просматривается на приведен- 
ных иллюстрациях, объясняется тем, что конфи- 
гурация изолиний отражает напряженно-дефор- 
мированное состояние, складывающееся из двух 
составляющих - статического, отвечающего пред- 
варительному поджатию образца, и динамическо- 
му, связанному с распространением возмущения 
от приложенного импульса. В связи с этим общий 
уровень напряжений в расчетной области растет, 
что на иллюстрациях смотрится как увеличение 
зоны возмущений в данный момент времени. 

2. Наличие в образце включения из другого мате- 
риала приводит к существенному изменению 
параметров напряженно-деформированного со- 
стояния в волне. На иллюстрациях с изолиния- 
ми отчетливо прослеживаются зоны больших 
градиентов интенсивности деформаций, совпа- 
дающие с границами раздела основного матери- 
ала образца и включения. Именно эти градиен- 
ты порождают смещения границ разделов ос- 
новного материала и включений, что, в резуль- 
тате, приведет к изменению амплитудных пара- 
метров электромагнитного отклика образца. 

Работа выполнена при поддержке гранта РФФИ робъа 
№ 05 - 01 - 98005 . 
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